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Abstract 

Multistable Levy motions are extensions of Levy motions where the stability index is allowed to 
vary in time. Several constructions of these processes have been introduced recently, based on 
Poisson and Ferguson-Klass-LePage series representations and on multistable measures. In this 
work, we prove a functional central limit theorem for the independent-increments multistable 
Levy motion, as well as of integrals with respect to these processes, using weighted sums 
of independent random variables. This allows us to construct continuous approximations of 
multistable Levy motions. In particular, we prove that multistable Levy motions are stochastic 
Holder continuous and strongly localisable. 
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1. Introduction 

Recall that a stochastic process {L{t),t > 0} is called (standard) a—stable Levy motion if 
the following three conditions hold: 

(Cl) L(0) = 0 almost surely; 

(C2) L has independent increments; 

(C3) L{t) — L{s) ~ Sa{{t — s)^/", /3, 0) for any 0 < s < f and for some 0<a<2, —1</3<1. 
Here S'Q,(cr,/5, 0) stands for a stable random variable with index of stability a, scale parameter 
a, skewness parameter f3 and shift parameter equal to 0. Recall that a governs the intensity of 
jumps. 

Such processes have stationary increments, and they are 1/a—self-similar, that is, for all 
c > 0, the processes {L(ct),t > 0} and {c^/"L(t),t > 0} have the same hnite-dimensional 
distributions. An a—stable Levy motion is symmetric when = 0. Stable Levy motions, and, 
more generally, stable processes have been the subject of intense activity in recent years, both on 
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the theoretical side (see, e.g. [l3|) and in applications [l2|. However, the stationary property 
of their increments restricts their use in some situations, and generalizations are needed for 
instance to model real-world phenomena such as financial records, epileptic episodes in EEG or 
internet traffic. A signihcant feature in these cases is that the “local intensity of jumps” varies 
with time t. A way to deal with such a variation is set up a class of processes whose stability 
index a is a function of t. More precisely, one aims at defining non-stationary increments 
processes which are, at each time t, “tangent” (in a certain sense explained below) to a stable 
process with stability index a(t). 

Formally, one says that a stochastic process {X{t),t E [0,1]} is multistable if, for almost 
all t E [0,1), X is localisable at t with tangent process X^ an a(f)—stable process. Recall 
that {X{t),t E [0,1]} is said to be h—localisable at t (cf. 0 , 8 ) , with h > 0 , if there exists a 
non-trivial process A/, called the tangent process of X at t, such that 


X{t + ru) -X{t) , 

hm- r- = X.[u), 

r\0 


( 1 ) 


where convergence is in finite dimensional distributions. 

Let D[0,1] be the set of cadlag functions on [0,1], that is functions which are continuous 
on the right and have left limits at all t E [0,1], endowed with the Skorohod metric ds 0]- 
X and have versions in i5[0,1] and convergence in ([T]) is in distribution with respect to d^, 
one says that X is h—strongly localisable at t with strong local form A{. 

In this work, we will be concerned with the simplest non-trivial multistable processes, 
namely multistable Levy motions (MsLM), which are non-stationary increments extensions 
of stable Levy motions. Two such extensions exist 


1. The field-based MsLM admit the following series representation: 

rF(i) = C'g'‘> l[„,„(X)Y<-'/”«» ((6|0,T|), 

(X,Y)Gn 


( 2 ) 


where H is a Poisson point process on [0,1] x M with mean measure the Lebesgue measure 
£, := sign(a)|a|^ and 


Cti = ( / X “sin(a;)da; 


-1 


( 3 ) 


Their joint characteristic function reads: 


f r+oo 

-2 / sin^ 

i[0,T] Jo 




{ m 

i''^^9jLF(tj) ^ = exp 

' ( 4 ) 

for d E N, ( 6 * 1 ,..., 9fi) E and {ti,... ,td) ^ These processes have correlated incre¬ 
ments, and they are localisable as soon as the function a is Hblder-continuous. 
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2. The independent-increments MsLM admit the following series representation: 


0(()= 5 ; ((6|0,T1), (5) 

(X,Y)Gn 


As their name indicates, they have independent increments, and their joint characteristic 
function reads: 


Eexp 


i 'y ^ 9jL]{tj 


= exp 


- lE 


^A[o, t 


](s) 


J=1 




a(s) 


ds 


( 6 ) 


for d G N, (01 ,..., 9d) G and (ti,... ,td) G These processes are localisable as soon 
as the function a verihes: 


(^a{x) — a{x + Inf —)■ 0 (7) 

uniformly for all x in hnite interval as t \ 0 j8|. 

Of course, when a{t) is a constant a for all t, both Lp and Lj are simply the Poisson representa- 



In this paper, we prove a functional central limit theorem for independent-increments MsLM: 
we show that certain weighted sums of independent random variables converge in (-D[0,1], ds) 
to Lp This allows us to obtain strong localisability of these processes. Moreover, we establish 
continuous approximations of MsLM and an alternative representation for the integrals of 
multistable Levy measure. Some properties of the integrals of multistable Levy measure are 
investigated. In particular, we prove that MsLM are stochastic Holder continuous and strongly 
localisable. 

The paper is organized as follows. In Section [21 we present the functional central limit 
theorem for independent-increments MsLM. In Section [21 we establish continuous approxima¬ 
tions of MsLM. In the last section, we give a representation of MsLM and investigate some 
properties, including stochastic Holder continuous and strongly localisable, of the integrals of 
multistable Levy measure. 


2. Functional Central Limit Theorems for Multistable Levy Motions 

We show in this section how to approximate the independent-increments MsLM in law by 
weighted sums of independent random variables. 

Theorem 2.1. Let {an{u))n,a{u),u G [0,1], he a class of cddldg functions ranging in [a,b] C 
(0,2] such that the seguence (a)„ tends to a in the uniform metric. Let (X(/c, n))^gpj 2 ^ 
he a family of independent and symmetric an{-^)—stable random variables with unit scale pa¬ 
rameter, i.e., X{k,n) (1,0,0). Then 
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• the sequence of processes 


Lf(u) = 


1 1 11 
1 \ 


E 

k=l 


X{k,n), 


U e [0,1], 


( 8 ) 


tends in distribution to Lj{u) in (-D[0,1], ds)-, where \x\ is the largest integer smaller than 
or equal to x. In particular, if a satisfies condition then Li{u) is localisable at all 
times. 


• the sequence of processes 


rp(u) 


[2^uJ 


k=l 


E 




X{k,n), 


u e [0,1], 


( 9 ) 


tends in distribution to Lj{u) in (-D[0, l],*^^), where (rj)j>i is a sequence of arrival times 
of a Poisson process with unit arrival rate and is independent of {X{k,n))^^^ 


• the sequence of processes 





X{k,n), 


u e [0,1], 


( 10 ) 


tends in distribution to Lj{u) in (Z1[0,1], 

Proof. We prove the first claim by the following three steps. 

First, we prove that L^p\u) converges to Li{u) in hnite dimensional distribution. For any 
Ui,U 2 G [0,1] and U 2 > Ui, we have, for any 6* G M, 


lim Ee' 

n—)-oo 




lim exp 

n—>-CX> 


L2"«2j 

fc=[2"tiij+l 


( 11 ) 


Notice that 

L2"n2j 


E 






L2"M2j 


< i'’ri°gi»i E ^ 


< 


A;=[2"uiJ+l 

aJ-) - a( 




a 


(12, 


where r = alp, i)(|6*|) + M[i,oo)(|6'|). By hypothesis, we have 


lim ||Q!n(-) - a(-)lloo = 0. 
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Thus inequality (IT^ implies that 


L2"«2j 

lim 

72^00 

k=[2'^ui\-\-l 




From fITT]) . it follows that 

lim 

n—>oo 


L2"«2j 

lim 

n—>-CX) 

k=l2^ui\+l 






(13) 


Hence L'f'\u 2 ) — L!f'\ui) converges in distribution and the characteristic function of its limit 
is defined by ([13]). Since L^P\u) has independent increments, the limit of L^P\u) has the joint 
characteristic function ([6|), i.e., L^P\u) converges to Lj{u) in finite dimensional distribution. 

Second, we prove that L^P\u) converges to Lj{u) in (Z1[0,1], ds). By Theorem 15.6 of 
Billingsley j^, it suffices to show that 


P 


Lf\u)-Lf\ui) > A, L^P\u2) 


W/ 




>AI<^ 

- 


U2 - Ml 


(14) 


for Ml < M < M2 iA > 0 and n > 1, where 7 = al[2,oo)(A) + 51(o,2)(A) and C is a constant 
depending only on a and b. If M2 — mi < 1/2", then either L^P\u 2 ) = L^P\u) or L^P\u) = 
L^P\ui)] in either of these cases the left side of f[T4l) vanished. Next, we consider the case of 
M2 — Ml > 1/2". Since L^P'\u) — L^p\ui) and L^P'\u 2 ) — L'f‘\u) are independent, it follows that 


p(|Li"^(M) - Lf\u^)\ > A, |4”^(m2) - Lf\u)\ > a) = 

p(|Lf)(M) - 4")(Ml)I > a) p(|l 4 (m 2) - Lf\u) 
Then, by the Billingsley inequality (cf. p. 47 of j^), it is easy to see that 


> A 


p(|l4(m)-l4(mi) 


> A < 


r 

^ A-2/A 

A 


(^1 - 


p 2IX / ( l2'^u \ 

/ ( 1 —p{- E 

d-2/A V L I 


■1^1 


«n(w) 


< 


'-2/A 

A 1 


'-2A 


E 


fc=[2"iiij+l 

\g\»n[^)d0 


fc=[2"uij+l 


de 
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C/A 

9 

/-2/A 



< y --1 9 'd 9 

- 2^2 
fc=[2"itij+l 

^ Cl r [_2^u\ - L2»mj 
“ [ 2^ 

where Ci is a constant depending only on a and b. Similarly, it holds 


P 


Lf(u,)-LT’{u) 


W/ 


> A < 


AT" 


[2”m2J - L2”mJ 


( 15 ) 


where C 2 is a constant depending only on a and b. Using the inequality xy < {x + y)‘^/4: for all 
x,y > 0, we deduce 

p(|l;”>(«) - L<”'(u.)| > A, |4 ">(u2) - 4”'(U)| > a) 


< 


< 


C 1 C 2 

A27 

^1^2 1 
4 A27 
1 

A^L 


[2'^mJ - [2"miJ 


[2”m2J - L2”mJ 


[2”n2j - L2 ”miJ 


1 2 


< C 1 C 2 


U 2 - Ui 


where the last line follows from the fact that 

[2”M2j - L2 ”miJ ^ 2^M2 - 2”mi + 1 


< 2 


U 2 - Ml 


This completes the proof of ffTTl) . 

Third, we prove that if a satishes condition 0. then Liiu) is localisable at all times. 
Falconer and Liu (cf. Theorem 2.7 of j^) have proved that the process Lj(u), dehned by the 
joint characteristic function ([H]), is localisable at u to Levy motions La{u){') with the stability 
index a{u). Here we give another proof to complete our argument. For any (ti, G [0,1]'^, 

from equality ([6]), it is easy to see that 


E exp <1 i 9j 
i=i 


Lj{u + rtj) - Lj{u) 


d 


= exp 


/ I «+rq](s) 

i=i 


a{s) 


ds > . 


Setting s = M + rt, we hnd that 

d 


E exp < i 'y^ 9j 


i=i 


Li{u + rtj) - Li{u) 

j,l/a(u) 
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= exp 



a{u-\-rt) 


^{a{u)-a{u+rt))/a{u) 


By condition it follows that 


linir^“^“^ "(“+^d)/“(“) — 1 \im. a{u + rt) = a{u). 

r\0 r\0 

Hence, using dominated convergence theorem, we have 

d 


( 16 ) 


lim E exp 

r\0 




9, 


Lj{u + rtj) — Lj{u) 


^l/a(u) 


= exp 


i=i 

d 

E exp < i E 9jL(y(^u) {t 


a{u) 


dt 


j) ( 1 


i=i 


which means that Li{u) is localisable at u to an a( m)— stable Levy motion La{u){t)- This com¬ 
pletes the proof of the first claim of the theorem. 

Next, we prove the second claim of the theorem. For any ui,U 2 G [0,1] and U 2 > ui, it is 
easy to see that, for any 6^ G M, 


n—>-oo 



(17) 


where the last line follows from the weak law of large numbers. Notice that (u) also has 
independent increments. The rest of the proof of the second claim is similar to the proof of the 
hrst one. For this reason, we shall not carry it out. 

In the sequel, we prove the third claim by the following two steps. 

First, we prove that L^^\u) converges to Lj{u) in hnite dimensional distribution. It is 
worth noting that does not have independent increments. This property implies that we 
cannot use the previous method. For any {ui,...,Ud) G [0,1]'^ and any (6*i,..., 6*^) G such 
that Q = uo < ui < U 2 < ■■■ < Ud, we have 


lim 

n—^oo 


d L2"Md d 

lim Eexp ^ i E E E«. 

1=1 j=l 

d L2"«iJ d 

\imEexp{-J2 E 

1=1 fc=[2"«i_ij j=l 


1 \ 


X{k, n) 


1 2" 


2 " To 


7 







= exp < - 


-E/ 

1=1 j=i 


|a(s) 


= exp 



ds 

a(s) 


^jrl[o,uj)('S) ds 
i=i 

which gives the joint characteristic function of Lj. 

Second, w^rove that L^^\u) converges to Lj{u) in (Zi)[0,1], ds). Again by Theorem 15.6 
of Billingsley |2| , it suffices to show that 

1 2 


P 


lP(u) 




>A, 


rc’('“2) 




>ai<A 

- A^t L 


U2 - Ml 


(18) 


for Ml < M < M 2 ,A > 0 and n > 1 , where 7 = al[ 2 ,oo)(A) + 61 (o, 2 )(A) and C is a constant 
depending only on a and b. We need only consider the case of M 2 — mi > 1 / 2 ”. Since L^^\u) — 
L^^\ui) and l[?^(m 2 ) — L^^\u) are conditionally independent given r 2 r., it follows that 


P 




. (n) 

"C 

77.) 


>A, 

p(|l7(«)-l7(“i) 




4 ” 7 ) 


> A 


r2n I — 


> A 


r 


2n 


> A 


It is easy to see that 


P 


4"'M-l7(«.) 


W/ 


> A 


2n J < 


A /•'/" 


A /•'/" 


'-2/A 


2 

■E 


p(|4”’(«2)-4‘’(«) 


r 


2n 


(19) 


-2/A 

exp 


E 


[2"tiJ 

E 

fc=[2"uij+l 


2" 




■1^1 




2" 


r 


2n 


dO 


< 


< 


A 1 


2 

Cl 2 


E 




271 


fc=[2"uij+l 

” r - [2 ”miJ 


x-yr 


2n 


where Ci is a constant depending only on a and b. Similarly, it holds 


P 


L7('“2) 


4 " 7 ) 



^ \ C2 2” r 

> A 

r 2 ™) < 

J - Xy T 2 n 


[2 ”m2 J - 


( 20 ) 


where C 2 is a constant depending only on a and b. From (IT^ . we find 


p(|l[?)(m) -L[?)(mi)| > A, |l[?)(m 2 ) -4"^(m) 


> A 


< 


< 


C 1 C 2 

■ [2”mJ - [2”miJ ■ 

■ [2”M2j - L2"mJ ■ 

E 

If Vl 

A27 

2n 

2” 


llraJ J 


C 

A^ 


n 2 


M2 — Ml 























































where C* is a constant depending only on a and b. This completes the proof of ([T5]) . 
Remark 2.1. Let us comment on Theorem \2.1\ 


□ 


1. We can define the independent-increments MsLM {Lj{x) : x E M.} on the whole line 
as follows. Let a{x), x E M, be a continuous function ranging in [a, 6 ] C (0,2], and 
satisfies condition uniformly for all x in finite interval as t \ 0. Set the functions 

= a{x + k) for all k > 0 and x E [0,1]. For any ak{x), by Theorem \2.1[ we can 
construct MsLM 

Ljfix) : [0,1] -)■ M, k>t). 

Taking a seguence of independent processes Lifix), x E [0,1], we define {Li{x) : a; > 0} 
by gluing together the parts, more precisely by 

yx\-i 

Li{x) = - [xj) + ^ -^ 4 ( 1 ), for all x > 0. (21) 

fc =0 

Similarly, for x < 0, we can define Lfix) = Lfi—x), since the function fi{x) = a{—x) is 
defined on [ 0 ,+oo). 

2 . Let (0(n))„gN be a seguence of numbers satisfying (pin) —)■ 00 as n —)■ 00 . Assume that 
a{u) is continuous in [0,1]. By an argument similar to the proof of Theorem \2.1[ the 
seguence of processes 



V«(^) 


X{k, n) 


tie [0,1], 


( 22 ) 


tends in distribution to Lj in (11[0, l],^^). Since a{u) is continuous, it is easy to see that 


a 


( [4>{n)u\ \ 

V </>(^) ) 


a{u) 


as n —)■ 00. 


By the fact that the summands of verify 


1 


X([ 0 (n)Mj,n) ~ S' 

<)>(n) > 



l/a( 


4>(n) } 



•} 


eguality means that the increment at the point u of an a{u) —multistable process 
Li{u) behaves locally like an a{u)—stable random variable, but with the stability index 
a{u) varying with u. 

3. If a{u) = a for a constant a E (0,2], then Li{u) is just the usual symmetric a—stable 
Levy motion Lain). Hence, ineguality ^2^) gives an eguivalent definition of the symmetric 
a—stable Levy motions: there is a seguence of independent and identically distributed 
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(i.i.d.) symmetric a—stable random variables (yfc)fcgN with an unit scale parameter such 
that 

\nu\ 

= (23) 

k=l ' 

tends in distribution to in (-D[0,1], ds)- This result is known as stable functional central 
limit theorem. 

4. A slightly different method to construct Lj{u) can be stated as follows. Assume that 
2 " *'5 a family of independent and symmetric a{^) —stable random vari¬ 
ables with the unit scale parameter. Then it holds 


[2"«J ^ X 

L,(u}= lim V (i-) 

^ ^ ^ \2^J \2^ 

k=l 


u e [0,1], 


(24) 


where convergence is in (-D[0,1], ds). To highlight the differences between the two methods 
fSf) and 1^]), note that X{-^) = while X{k,n) and X{2k,n-\- 1) are two i.i.d. 

random variables. 

5. Inspecting the construction of field based MsLM in Falconer and Levy Vehel Q/, it seems 
that the seguence of processes 


lP{u)= (^) ' ^ e (25) 

k=\ 

tends in distribution to Lf{u) in {D[0,l],ds). Unfortunately, it is not true in general. 
We have the following counter example. 

Example 1. Consider the case of an{u) = a{u) = |l{o<u<i} + charac¬ 
teristic function of L^fi\u) is given by the following eguality: for any 6^ G M, 


[2"mJ 




n 


k=l 


_ ^ |0|«(2.)(_) { mg [0,1], 


Since, for all u G (|, 1] and 6^0, 


[2"nJ 




^ / 1 \ «(^)/"(“) 


[2"fe/2j 


> 


k=l 


k=l 


\ b/2u 


(26) 


oo, n ^ oo, 


(27) 


we have L^fi\u) —>■ 0 for all u G (|,1]. 
Lf{u) in (Z1[0, l],(i 5 )- 


Thus 


lP{u) 


does not tend in distribution to 
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3. Continuous Approximation of MsLM 


It is easy to see that when a{u) is a constant, then the independent-increments MsLM 
reduce to a—stable Levy motions. It is well known that a—stable Levy motions are stochastic 
Holder continuous but not continuous. We wonder if there exists a continuous approximation 
of independent increments MsLM? The answer is yes. 

3.1. A continuous stable process 

First, we shall construct a continuous stable process. To this end, we shall make use of the 
following useful theorem. 

Theorem 3.1. If the i.i.d. random variables {Zjk)j^k follow an a—stable law, then it holds, for 
all c> 1 /a, 



P 


Proof. We only need to show that, for all c > 1/a, 



By Borel-Cantelli Lemma, it is sufficient to prove, for all c > 1/a, 

^^P f max \Zjk\ > < oo. (28) 

j>l V^=0v-,2J-l J 

To prove fl2^. we need the following technical lemma (cf. Property 1.2.15 of Samorodnitsky 
and Taqqu [l3[ for details). 

Lemma 3.1. Let Z ~ Sa{(T,ld,p) with 0 < a < 2. Then 


hm;,^o,A“P(Z> A) = 


limA^oo A“P(Z <-A) = 


Return to the proof of fl28l) . For all c > 1/a and all j large enough, we have 



max 



1 - 


n I* i r°) ■ 


(29) 


fc =0 


11 


Then, by equality (12^ and Lemma 13.11 we deduce 


P ( max \Z~k 
\k=0,...,2j-l ^ 



1- (i + o(^)) 

O ( 7-7-—TT^ , ?' —)■ CXO 


Thus we obtain fl28l) for all c > 1/a. □ 

In the following theorem, we give a construction of continuous stable process. First, we 
recall the definition of the “triangle” function: 

( 2t forte [0,1/2) 

^{t) = <^ 2 - 2t for t e [1/2,1] 
y 0 otherwise. 

Dehne <^jfc(t) = ^p{2H — k), for j = 0,1,..., and fc = 0,..., 2-^ — 1. 

Theorem 3.2. Assume the i.i.d. random variables follow a symmetric a—stable law 

with the unit scale parameter. Then, for all d> 1/a, the process 

00 2^-1 

= E E t e [0,1], 

j=0 k=0 

is a continuous and symmetric a-stable process. When d = 1/a, the process X(t) is also a 
symmetric, may not be continuous, a-stable process in U’{VL x [0,1]) for any D < p < a. 

Proof. Set X_i = 0 and define the sequence of processes {Xj)j^^ by: 

2J-1 

~ 2 ^ Zjk^pjkit). 

k=0 

First we show that the sequence of processes {Xj)j^fq converges almost surely uniformly. Indeed, 
for all t, 

2^-1 

X^{t) - X,_,{t) = 2-^^Z^kTjk{t). 

k=0 

Since the functions {(pjk)j,k have disjoint supports and \(pjk\ < 1, it follows that 

\\X,{t)-X,_,{t)\\^ = 2-^'" max \Zjk\. 

k=0,...,2^ ^ 

Theorem 13.11 entails that (Xj)jgN converges almost surely in (^([O,!], || • ||oo) to a continuous 
process X for all d > 1/a. When d = 1/a, we show that the sequence (Xj)jgN converges to a 
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random variable X in L^{Q x [0,1]) for any 0 < p < a. Indeed, for any 0 < p < a, 

pi 2^-1 pi 

Jo Jo 

< 2-^^>/“E|Zoor [ vloit)dt 
Jo 

= 2-^'P/“E|Zoo|^, 


( 30 ) 


this entitles convergence of {Xj)j^^. 

Next, we prove that X is a symmetric a—stable process. By Theorem 3.1.2 of Samorodnitsky 
and Taqqu (1994), we only need to check that all linear combinations 

d 

'^bkXitk), d>l, fi, ...,td e [0,1] and & 1 , ...,6rf real 

k=l 

are symmetric a—stable. We distinguish two cases as follows. Define 

Dr, = I —:0<k <2^ 


and D = Un=o,i,...^n- 

i) If tk € D, then all random variables X{tk),l < k < d, are symmetric a—stable. Thus all 
linear combinations Ylt=i^kX(tk) are symmetric and a—stable. 

ii) For tk E [0,1], 1 < /c < d, we have tki G D such that tki tk,l ^ oo. Since X is continuous, 
we have 


'^bkXitk) = lim '^bkXitki). 

k=l k=l 

Its characteristic function has the following form: 

d \ p' d 


Eexp li6^bkX{tk)\ = hm Eexp bkX{tki] 

^ k=l ^ ^ k=l 


1/a 


It is easy to see that the scale parameter of X]fc=i ^kJ^dki) is 

( oo 2i-l d 

|6fc|2"^V(2^4z - i) 

j=0 i=0 k=l 

Since at most one summand of the sum 2~^'^ip{2H — i) is non-zero and 

2^-1 d 

Y, (5^ |6fc|2-^V(2'4z - *))“ < d6“2-^“^ 

i=0 k=l 
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where b = max{\bk\, 1 < k < d}, then a{a) = lim^^ooexists for d > 1/a and 
Eexp I iO bkX (t^) | = hm exp | — cr;(a)"|6'|“| 

= exp I — (T(a)“|6'|"|. 

This implies that all linear combinations 'Yfk=ibkX{tk) are symmetric a—stable random vari¬ 
ables. This completes the proof. □ 

One deduces the scale parameter a{t) of the process X{t) is given as follows 

oo 2>-l 

-“w = EE( 2-^‘^ip{2H - k)y. 

j=0 k=0 

By noting that at most one ip{2H — k) is non-zero for all j, we have the following estimation of 
the scale parameter 

//“(t) < a{t) < ’ ^e[ 0 ,l]. 

It is worth noting that when t 7 ^ 0,1, we have a{t) > 0. This observation will be useful to 
establish continuous approximations of MsLM in the next subsection. 


3.2. Continuous approximations of MsLM 

In Theorem 12. we establish discrete approximations of the independent-increments MsLM. 
In this subsection, we shall give continuous approximations of the independent-increments 
MsLM. It is worth to noting that one cannot make use of the method of Theorem 13.21 to estab¬ 
lish continuous approximations of MsLM in general, since a sum of two stable random variables 
with different stability indices is not a stable random variable. To obtaining continuous approx¬ 
imations of the independent-increments MsLM, our main method is to replace the summands 
in (jH]) by a sequence of independent and continuous stable processes starting at 0 , for instance 
the stable processes established in Theorem 13.21 

Theorem 3.3. Let a{u) be a continuous function ranging in [a, 6] C (0,2]. Assume that 

( X,j^Af) ) is a family of independent and continuous a(X)—stable random pro- 

V VnGN,fc=0,...,2"-l ^ 

cesses. Assume = 0 and cr^(^_^^{t) > 0 for all t e (0,1] and alln eN, k = 0,..., 2 "' — 1 , 

where is the scale parameter of X^^_k_s^{t). Define 
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Then (S'„)„gp} is a sequence of continuous processes and the process Sn{u),u G [0,1], tends in 
distribution to Lj{u) in (D[0,1], (i^). 

By the definition of 5 '„(m) in (1^ . it seems that the process Sn{u) restores more and more 
details of Li{u) when n is increasing. 

It is worth noting that when a{u) = a for a constant a G (0, 2], Theorem l3.3l gives continuous 
approximations to the usual symmetric a—stable Levy motion La{u). 

Proof. It is easy to see that the first item in the right hand side of flMl) converges to zero in 
distribution as tt, —)■ oo, i.e., 


lim 

n^oo 


( 1 





in distribution. Notice that the summands 


a(4r)(2") 


a, 


I - 


1 

2n 


0 


(32) 


(33) 


in the right hand side of fl31l) are independent a(^)—stable random variables with the unit scale 
parameter. Using Theorem 12.1[ we find that the process Sn{u),u G [0,1], tends in distribution 
to Lj{u) in (T)[0,1], ds). □ 


4. Integrals of Multistable Levy Measure 

Let a = a{u),u G [0,1], be a cadlag function ranging in [a, h] C (0, 2]. Denote by 

£q,[0 , 1] = |/ : / is measurable with ||/||a < oo|, 


where 


Note that 


:= inf < A > 0 : 




A 


a.{x) 


dx = 1 


and 


= 0 . 


\a is a quasinorm; see Falconer and Liu and Ayache [i|. Using the Kolmogorov 
consistency conditions and the Levy continuity theorem, Falconer and Liu [sl (see also Falconer 
(il) proved that the characteristic function, for all (di, ...,6^) G M'^, 


Eexp I* ( ^ I“ / I S 


a{x) 


dx 


(34) 


well defines a consistent probability distribution of the random vector (/(/i), L(/ 2 ), •••, Hfd)) ^ 
on the functions fj G £q[ 0, 1], where /(/) = J f{x)Ma{dx). They called Mq, the multistable 
Levy measure and I{f) = J f{x)Ma{dx) the integral with respect to Ma- Moreover, they also 
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showed that the integrals of functions with disjoint supports are independent. In particular, it 
holds 

Li{u) = j l[o, „](x)M„(dx), u e [0,1]. 

In the following theorem, we give an alternative dehnition of the integrals based on the 
weighted sums of independent random variables. 

Theorem 4.1. Let {X{k,n))^^^ 2 ^ defined by Theorem \2.1{ Then, for any f G 

£q,[0 , 1], it holds 

^ / 1 \i/«(^) / k\ 

/ Et) (35) 

k=i 

in distribution. 


Proof. Denote by 

k 2 ^ 

^ arid Xn = ^S{k,n). 

k=l 

It is easy to see that, for any 6 * G M, 

Ee^^^- = = exp 

k=l 



Hence, we have 


lim 

n^oo 


exp 






which means lim„_j.oo^n = Jq fix)Ma{dx) in distribution by the definition of the multistable 
integrals with respect to the multistable Levy measure M„. □ 

The following theorem relates the convergence of a sequence of a( m)— multistable integrals 
to the convergence of the sequence of integrands. 


Theorem 4.2. Assume Xj = fj{x)Ma{dx) and X = f{x)Ma{dx), for fj, j 
£q,[0,1]. Then 

lim Xj = X 

j^oo 


in probability, or 
in distribution, if and only if 


lim {Xj-X) = 0 

j^oo 


lim Wfj- fWa = 0. 


3^00 


1,2,...J 


G 
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Proof. The convergence linij_^oo ^ in probability is equivalent to linij^oo(-^i ~ ^) = 0 in 
probability and hence to the convergence in distribution to zero of the sequence {Xj —X)j=i^ 2 ,...- 
If Xj — X convergence in distribution to 0, then, for any 6^ G M, 


1 = lim = lini 


J-J-CXD 


3^00 


exp 


- fix) 


a(x) 


dx 


(36) 


which is equivalent to, for any A > 0, 


lim 

j^oo 



fj{x)- f{x) 

lo 

A 


dx 


0 . 


This equality means linij^oo ll/j-/ll« = 0. □ 

The last theorem shows that convergence in probability of multistable integrals coincides 
with convergence in quasinorm || ■ ||q,. 

The convergence hmj_j.oo Xj = X almost surely implies the convergence limj_j.oo Xj = X in 
probability. Thus the following corollary is obvious. 


Corollary 4.1. Assume that Xj, j 


1,2,.. and X are defined by Theorem 4.2 


If 


lim Xj = X 

j^QO 

almost surely, then 

lim Wfj- fWa = 0. 

4 . 1 . Independence 

Independence of two multistable integrals imposes a stronger restriction on the integrands: 
they must almost surely have disjoint supports with respect to Lebesgue measure C. Indeed, 

Theorem 4.3. Let Xi = fi{x)Maidx) and X 2 = f 2 ix)Ma{dx) be two multistable inte¬ 
grals, where fj G £o[0,1], j = 1,2. Assume either 

[a, h] C (0, 2) 


or 


fi{x)f 2 {x) > 0 C — a.s. on [0, 1]. (37) 

Then Xi and X 2 are independent if and only if 

fi{x)f 2 {x) = 0 C — a.s. on [0, 1]. (38) 
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Proof. Two multistable integrals Xi and X 2 are independent if and only if, for any ( 6 ^ 1 , 62 ) G 


Eexp |i( 6 *iXi + 6 * 2 X 2 )! = Eexp |i 6 *iXi| Eexp |i 6 * 2 X 2 


( 39 ) 


Notice that 


Eexp i( 6 *iXi + 62 X 2 ) \ = exp 



1 , 2 


and that 


t=i 

2 rl 


X 


a{x) 


dx > , 


Eexp |i 6 *iXi| Eexp |* 6 * 2 X 2 ! = exp J 9jfj{x) 


a{x) 


dx 


Equating the moduli of fl5^ gives 

n , 2 


^Ojfj{x) ^^dx = J 2 [ 

j=i j=i Jo 


a{x) 


dx. 


(40) 


Notice that fl40|) implies that 


fl{x) - f 2 {x) 


a{x) 


dx = 




a{x) 


dx + 


f 2 {x) 


a{x) 


dx. 


fl{x) + f 2 {x) 


a(x) 


dx 


(41) 

(42) 


Assume [a, b] C (0, 2). We argue as Lemma 2.7.14 of Samorodnitsky and Taqqu [l^. When 
a G (0, 2), the function ra{u) = n > 0, is strictly concave. Therefore, for £x x G [0,1], 

|/l(x)+/2(x)|“(") + |/l(x)-/2(x)r(") 

raix)(,\fi{x) +/2((r)n +r„(^)(|/i(x) -/2(x)n 


= 2 


^ 2 


1/1 (a;) + f2{x)\^ + 1/1 (a;) - /2(a()p 


= 2 r„(,)(/i(i)^ + / 2 (i)^) 

< 2 (i/,(i)rw + i/ 2 (ni“'">) 

with equality in the preceding relations equivalent /i(x)/ 2 (x) = 0. Inequalities fITTD and 
imply that 


(43) 


h(x) - hix) 


a(x) 


dx + 


fl{x) + f 2 {x) 


a{x) 


= 2 


fiix] 

18 


a(a:) 


dx + 


dx 


f2{x) 


a{x) 


dx 1 . 


(44) 


























Now implies that the left-hand side of the last equality is always less than or equal to the 
right-hand side of the last inequality and, if they are equal, then necessarily fl55]) holds. 

Assume ([37]). Then it holds \fi{x) — f 2 {x)\ < \fi{x) + f 2 {x)\ C — a.s. on [0, 1]. When 
a G (0, 2], the function ra{u) = is increasing in u G [0, cxo). Hence flT2l) holds if and only if 
fl38|) holds. 

This proves that fl55]) is a necessary condition for the independence of Xi and X 2 . It is also 
sufficient because if fl55]) holds, then (HUD also holds. □ 

The preceding result is very useful and will often be used in the sequel. 

Theorem 4.4. Assume fj G £q-[0, 1], j = 1,..., d. Assume either [a, h] C (0, 2) or fi{x)fk{x) > 

0 C—a.s. on [0, 1] for any subset {i, k} o/{l, 2,..., d}. The multistable integrals Xj = fj{x)Ma{dx), 
j = 1,..., d, are independent if and only if they are pairwise independent, i.e., if and only if 

fi{x)fkix) =0 £ - a.s. on [0, 1] (45) 

for any subset {i, k} of { 1 , 2 ,..., d}. 

Proof. Independence clearly implies pairwise independence. By Theorem 14.31 pairwise inde¬ 
pendence implies fl45D . If (145]) holds, then it holds, for any {9i,...,9d) G M'^, 



Thus the joint characteristic function of Xi, ...,Xd factorizes 



d 

= Eexp f^i9jXj I. 


This proves that Xi, ...,Xd are independent. □ 

4 . 2 . Stochastic Holder continuity 

We call a random process X{u), u E I, is stochastic Holder continuous of exponent (3 G (0,1], 
if it holds 

limsup P(|X(m) — X{r)\ > C\u — r\^) = 0 

u^rGl, |ti—r|^0 

for a positive constant C. It is obvious that if X{u) is stochastic Holder continuous of exponent 
(3i G (0,1], then X{u) is stochastic Holder continuous of exponent (32 G (0,/5i]. 
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Example 2. Assume that a random process X{u), u & I, satisfies the following condition: 
there exist three strictly positive constants 7 , c, p such that 


E\X{u) — Xir)^ < c\u — rf, u,r E I. 


Then X{u),u G I, is stochastic Holder continuous of exponent ft G (0, min{l, p/ 7 }). Indeed, it 
is easy to see that for all u,r E I, 


f(^\X{u) - X{r)\>C\u-r\^'^ < 

< 


E|X(m) -X(r)|T' 
C'y\u - r\h'y 


c 


u — r 


P-0'1 

5 


which implies our claim. 


The following theorem gives a sufficient condition such that the integrals with respect to 
multistable Levy measure are stochastic Holder continuous. 

Theorem 4.5. Assume that X{t) = f{t,x)Ma{dx) is a multistable integral, where f{t,x) 
is jointly measurable and f{t,x) E £q,[0, 1] for all t E I. If there exist two constants rj > 0 and 
C > 0 such that 


Then it holds 




a{s) 


L 

f{'t,s) - f{v,s) 

ds < C 

t — V 


(47) 


V{\X(t) - X(t,)| > |i - t,|f) < \t - t,v£l, 


(48) 


where Ca^ is a constant depending on a, b and C. In particular, it implies that X(t) is stochastic 
Holder continuous of exponent fl E (0, min{l, 77 / 6 }). 

Proof. By the Billingsley inequality (cf. p. 47 of j^) and fl47|) . it is easy to see that, for all 
t,v E I and all a: > 0, 


P(|X(f) -X(n)| > x) 
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Taking x = \t — np, we obtain fl48p . This implies that X{t) is stochastic Holder continuous of 
exponent (5 G (0, min{l, 77 / 6 }). □ 

As an example to illustrate Theorem 14.51 consider the weighted MsLM introduced by 
Falconer and Liu j^. The following theorem shows that the weighted MsLM are Holder con¬ 
tinuous of exponent f3 G (0, min{l, !/&}). 

Theorem 4.6. Let ^ 

Y{t)=[ w{x)lio^ t]{x)Ma{dx), tG[0,l], 

Jo 

be a weighted multistable Levy motion, where the function w{x),x G [0,1], is cddldg. Then Y(t) 
is stochastic Holder continuous of exponent (J G (0,min{l, 1/b}). Moreover, it holds 

P(|y(t)-F(t;)| > \t-vf)<Ca,b [0,1], (49) 

where Ca,b is a constant depending on a, h, «(•) and w{-). In particular, it implies that Lj{u),u G 
[0,1], is stochastic Holder continuous of exponent f3 G (0,min{l, 1/&}). 

Proof. Set f(t, x) = ia(a:)l[o,i](x), t, x G [0,1]. It is easy to see that, for all v,t E [0,1] such that 

V < t, 




a{s) P 



/(^,s) - f{v,s) 

ds< 

w;(s)l[„, t](s) 

Jo 


Jo 

P 



VI 

i 1 [d, t] ("S) ds 
'0 


< Cujit-v), 


where = sup^gjQ j^] By Theorem 14.51 we get fl49]) . This completes the proof of 

Theorem 14.61 □ 


4-3. Strongly localisability 

When the function a{x) G [a,b],x E [0,1], is continuous, some sufficient conditions such that 
the multistable integrals are localisable (or strongly localisable) has been obtained by Falconer 
and Liu. In the following theorem, we give some new conditions such that localisability can be 
strengthened to strongly localisability. 


Theorem 4.7. Assume that f{t, x) and h(t, x) are jointly measurable; and that f(t, x), h{t, x) G 
£q,[0, 1] for any t E [0,1]. Assume that X{t) = f(t, x)Ma{dx) and X'^{t) = h(t, x)Ma(dx) 
are two multistable integrals and have versions inD[0, 1]. Suppose thatX(t) is 1 /a{x) —localisable 
at X with local form X'^{t). If there exist two constants rj > 1 and C > 0 such that 


f{x + rt, s) — f{x + rv, s) 


..l/a{x) 


a{s) 


ds < C 

t — V 


t,v e |o, 1], 


(50) 


for all sufficiently small r > 0, then X{t) is strongly localisable at all x E [0,1]. Moreover, if 
X{f) has independent increments and fl50|) holds for a constant 7] > 1/2, then the claim holds 
also. 
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Notice that condition fISUI) is slightly more general than the condition of Falconer and Liu 
(cf. Theorem 3.2 of jsl): there exist two constants rj > 1/a and C > 0 such that 


f{x + rt, •) — f[x + rn, •) 


rf.l/a{x) 

for all sufficiently small r > 0. 
Proof. For any x G [0, 1), define 

Xr{u) = - 


< c 

t — V 

OL 


-X{x) 



t,v € |0,1], 


(51) 


j.l/a{x) 


By Theorem 15.6 of Billingsley j^, it suffices to show that, for some /3 > 1 and r > 0, 


P 


Xr{u) — Xr{Ui) > A, Xr{u2) — Xr{u) 




U 2 — Ml 


(52) 


for Ml < M < M 2 , A > 0 and r G (0,1], where (F is a positive constant. Since ^^(m) — ^^(mi) and 
Xj.{u 2 ) — Xr{u) are symmetric, it follows that 


P 


>A, 


X.f{xi2) X.f.{xPj 


> A 


< 


Xr{u) — Xr{Ui) 

4p(^X,(m) -X,(mi) + {Xr{u2) - Xr{u)^ > 2A 

= 4P^Xj.(m 2) — Xr(Mi) > 2Aj. 

By the Billingsley inequality (cf. p. 47 of j^) and (15(11) . we have 

p(x..(ti2) - X..(tii) > 2 a) - 

J-l/X 


p1/X pI 

< X 

J-l/X Jo 

nl/X 

= A / W 


f{x + rM 2 , s) - f{x + rMi, s) 


a{s) 


r^/a{x) 

^ f{x + ru 2 ,s)-f{x + rui,s) 


< — 
- X-f L 


J-l/X Jo 

Cl r I*? 


r 


l/a{x) 


ds dO 

a{s) 


ds d6 


U2 — Ml 


(53) 


where fi = al[i^oo)(6') + &l(o,i)(^),7 = al[i,oo)(A) + 61(o,i)(A) and Ci is a positive constant 
depending only on a, b and C. Thus 


P 


Xr{u) — Xr{Ui) > A, Xr{u2) — Xj,(m) 


> A < 


4^1 

A"^ 


M2 — Ml 
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Hence, by (IS2]) . if r; > 1, then X(t) is 1/a(x)—strongly localisable at x with strong local form 
If X [t) has independent increments, then 


P 


Xriu) — XriUi) 


> A, 


Xr(^^2} Xj<{Vj^ 


> A 


= P 


Xr{u) — Xr{Ui) 


> A P 


Xr{u 2 ) - Xr{u) 


> A 


(54) 


By an argument similar to (I53p . it follows that 


P 


Xriu) — XriUi) 


> A < 


4Ci 

X I 


u — Ui 


and 


P 


Xriu^^ Xriu) 


>A < 


4Ci 

A-^ L 


U2—U 


Using the inequality xy < ix + y)'^/4, x,y > 0, we have 


P 


Xriu) — XriUi) > A, Xriu2) — Xriu) 


> A < 


16Ct 

A27 


2 r 


U 2 — Ui 


2 ri 


(55) 


Thus, if 2?7 > 1, by (l52|) . then X(t) is 1/a(x)—strongly localisable at x. This completes the 
proof of the theorem. □ 

As an example to illustrate Theorem 14.71 consider the weighted MsLM. Falconer and Liu 
have proved that the weighted MsLM are localisable. The following theorem shows that the 
weighted MsLM are not only localisable but also strongly localisable. In particular, it shows 
that the independent-increments MsLM are strongly localisable. 


Theorem 4.8. Assume that the function aiu),u G [0,1], satisfies condition Let 

Yit)= f wix)lio, t]ix)Maidx), tG[0,l], 

Jo 

he a weighted multistable Levy motion, where the function wix),x G [0,1], is continuous. Then 
Yit) is 1/aix)—strongly localisable at all x G [0,1] with strong local form wix)La{x)i-)■ In 
particular, this implies that L/(t) is 1/aix) —strongly localisable at all x G [0,1] with strong 
local form La(^x)i-), o,n aix) —stable Levy motion. 

Proof. It is known that Yit) is 1/a(x)—localisable at all x with strong local form wix)La{x)i-)] 
see Falconer and Liu j8|. Set fit,x) = tc(x)l[o, t]ix), t,x E [0,1]. By flTB]) . the integrand of 
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Y(t) satisfies, for all t,v G [0,1] such that v <t, 



f{x + rt, s) — f{x + rv, s) 

a{s) 

ds = 

Jo 

('^) 1 [x-\-rv, x-\-rt] ('^) 

L 

j,l/aix) 

j,l/a(x) 


a{s) 


ds 


w{x + rz)l[y^ t]{z) 


OL{x-\-rz) 


^(a(x)-a{x+rz))/a(x) 


< a, 




< 2 C^ {t - v), 

for all sufficiently small r > 0, where s = x + rz and = sup^gjo.i] \w{z)\°‘^^\ By the fact that 
the integrals of functions with disjoint supports are independent, it is easy to see that Y{t) has 
independent increments, the first claim of the theorem follows by Theorem 14.71 In particular, 
since Li{t) = l[o, t](x)-^a(dx), the first claim of the theorem implies the second one with 
w(x) = l,x e [ 0 , 1 ]. □ 


Remark 4.1. By inspecting the proof of Falconer and Liu l3j, we can see that Y{t) is also 
1/a{x) — localisable at allx with strong local form w{x)La(x)it) when the function w{x),x G [0,1], 
is cddldg. Hence, Theorem\4.8\ holds true when the function w{x),x G [0,1], is cddldg. 


References 

[1] Ayache, A. (2013). Sharp estimates on the tail behavior of a multistable distribution. 
Statist. Probab. Letter 83, 680-688. 

[2] Billingsley, P. Convergence of Probability Measures. John Wiley: New York, 1968. 

[3] Falconer, K. J. (2002). Tangent helds and the local structure of random helds. J. Theoret. 
Probab. 15, 731-750. 

[4] Falconer, K. J. (2003). The local structure of random processes. J. London Math. Soc. 
67(3), 657-672. 

[5] Falconer, K. J. (2012). Localisable, multifractional and multistable processes. Seminaires 
and Congres 28, 1-12. 

[6] Falconer, K. J., Le Guevel, R., Levy Vehel, J. (2009). Localizable moving average sym¬ 
metric stable and multistable processes. Stochastic Models, 25(4): 648-672. 

[7] Falconer, K. J. and Levy Vehel, J. (2009). Multifractional, multistable, and other processes 
with prescribed local form. J. Theoret. Probab. 22: 375-401. 


24 













[8] Falconer, K. J. and Liu, L. (2012). Multistable Processes and Localisability. Stochastic 
Models 28: 503-526. 

[9] Le Guevel, R., Levy Vehel, J. (2012). A Ferguson-Klass-LePage series representation of 
multistable multifractional processes and related processes. Bernoulli 18, No. 4, 1099-1127. 

[10] Le Guevel, R., Levy Vehel, J. (2013). Incremental moments and Holder exponents of 
multifractional multistable processes. ESAIM: Probab. Statist. 17, 135-178. 

[11] Le Guevel, R., Levy Vehel, J., Liu, L. (2012). On two multistable extensions of sta¬ 
ble Levy motion and their semimartingale representations. J. Theoret. Probab., DOI 
10.1007/sl0959-013-0528-6. 

[12] Nolan, J., Bibliography on stable distributions, processes and related topics, 
http://academic2. american.edu/jpnolan/stable/StableBibliography.pdf, 

[13] Samorodnitsky, G. and Taqqu, M. S. (1994). Stable Non-Gaussian Random Processes, 
Chapman and Hall, London. 


25 


